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ABSTRACT 

Problems where impact broadened radio recombination lines appeared nar- 
rower than predicted first showed up ~ 40 years ago at frequencies below ~ 3 
GHz. But it was soon found that the observations could be explained by throw- 
ing out the uniform density models and replacing them with variable density 
ones. However, this problem re-appeared recently when a mysterious line nar- 
rowing above quantum numbers of (n, An) = (202,8) was reported from sensitive 
observations of Orion and W51 near 6 GHz. Here it is demonstrated that the 
narrowing is unlikely to be caused by the data processing technique and therefore 
must be source related. It is further demonstrated that the observed line nar- 
rowing can be tied to one of the fundamental properties of radio recombination 
lines; namely the fact that the spacing of adjacent ri-transitions increases with 
frequency. The line narrowing is observed to begin when the n-transition den- 
sity, Dn, exceeds ~ 11.6 transitions per GHz. This may imply that it is somehow 
related either to a previously overlooked effect in the impact broadening process, 
or to some unknown parallel process, that is tied to the separation between ad- 
jacent n-transitions. Based on these results it can be concluded, as has also been 
concluded in several theoretical investigations, that the observed line narrowing 
is not tied to a fixed range of either n or An. 

Subject headings: HII regions~ISM: abundances-radio lines: ISM 



1. Introduction 

Using radio recombination lines (RRLs) 
we can obtain much new information 
extending all the way from the physi- 
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cal conditions inside galaxies, down to 
the details of processes that occur in- 
side individual atoms. As was discussed 



by iGordonI (120081 ) they can thus be a 
valuable tool for astronomers and physi- 
cists. That it would be possibl e to detect 
RRLs was first suggested by iKardashev 
( 119591 ). The theory of impact broa dening 
of the radio lines was developed by iGriem 



1 



(Il967l ). A broade r coverage of RRLs was 



carrie d out by iGordon and Sorochenko 



( 20021). However, it soon becarn e apparent 
f Brocklehurst and Seaton 1972h that there 



were problems. iHjellming and Davies 



(jl970l ) had found for several nebulae that 
for homogeneous models, an electron den- 
sity between A^e ~ 2.5 x 10^ and 1 x 10^ 
is required to explain the obse rved recom- 
bination line intensities. But Churchwell 



(jl97ll ) found much less broadening be- 
low ~ 3 Ghz than was predicted for this 
density and concluded that the widths 
predicted by Griem's theory appeared 
to be over-estimated at these low fre- 
quencies. It was concluded that for uni- 
form density models the observational 
results for the line-to-continuum ratios 
were not compatible with the l ine profiles 
(iBrocklehurst and Seatodll972l ). However, 
it was argued that the failure to detect 
broadening at low frequencies {u < 1.6 
GHz), and the smaller amount of ex- 
cess broadening detected at < 5 GHz 
than expected from a simple application 
of Griem's formula, were explicable in 
terms of density fluctuations in Orion 
( IChurchwelll Il97ll ). and models in which 
the electron density varied with radius 
were devised to explain the observations. 

2. More Recent Results 



2.1. Data at 6 GHz 

The low-frequency result discussed above 
was not the only place where mysterious, 
low impact broadening results were ob- 
tained. Using a novel, multiple-overlap 
frequency-switching observing technique 
we were able to observe for the first time 



of Rydberg-Rydberg transitions in Orion 
and W51 near 6 GHz that ranged from 
(n,An) = (102,1) to (289,25) and were al- 
most continuous in increasing An. When 
observations are carried out in a relatively 
narrow window and therefore close to the 
same frequency, the transitions available 
must all have different An values. The 
(n, An) transitions included in our win- 
dow were (102,1), (129,2), (147,3), (174,5), 
(184,6), (194,7), (202,8), (210,9), (217,10), 
(224,11), (230,12), (236,13), (241,14), (247,15), 
(252,16), (257,17), (261,18), (266,19), (270,20), 
(274,21), (278,22), (282,23), (286,24), and 
(289,25). 

Because we used the frequency switch- 
ing observing technique to obtain fiat base- 
lines without the need to fit and remove 
sinusoids or polynomials, the linewidths 
were affected by the data processing and a 
special procedure ha d to be used to recover 



the t rue hnewidths (IBelll 119971 iBell et al. 



20001). To obtain the true linewidths we 



first generated test Voigt profiles based 
on accepted impact broadening theory for 
each An-value, using the known parame- 
ters of each source, such as temperature 
and electron density. We then processed 
these test lines in a manner identical to 
the way the line profiles obtained at the 
telescope were processed. These processed 
test linewidths were then compared to the 
linewidths obtained for the processed tele- 
scope lines. As further evidence that our 
technique is likely to be as reliable as 
position-switching for strong sources like 
Orion it is noted here that the (Hel71,7) 
line was detected in Orion at the correct 



flBelll 119971 : iBell et al.ll2000L l201lf ) a series 



level of 3mk by iBelll ( 120001 ) using frequency 
switching, b u t was not detected in Fig 4c of 



Bania et al.l (Il987( ) using position switch- 
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ing. These authors attribute this failure to 
large systematic effects present in the mea- 
surement of lines in this strong continuum 
source when position switching is used. 

Observed line widths obtained for Orion 
and W51 at 6 GHz after processing are 
plotted vs n iri Fig 1. Data are from 
Bell et al.l ( 1201 ll ). The curve in Fig 1 rep- 
resents the linewidths obtained after pro- 
cessing the test lines. Below n = 200 
the processed telescope lines agree with 
the processed test lines indicating that the 
linewidths calculated using the impact the- 
ory give a reasonable fit to those observed 
at the telescope. Because the widths of 
the lines are significantly reduced by our 
reduction process it is important to note 
that the actual linewidths are much larger 
as c an be seen f rom t he upper curve in Fig 
4 of iBell et al.l fl2000l ). 

Fig 1 thus shows that the widths of the 
test lines and the telescope lines increase 
similarly with {n,An) up to (202,8), and 
exactly as predicted from theory. Above 
(202,8) the processed test lines showed no 
sign of decreasing in width. However, the 
widths of the lines obtained at the tele- 
scope can be seen to decrease abruptly for 
both the Orion and W51 data. The rel- 
evant result here is the fact that the tele- 
scope profiles narrowed abruptly above n = 
200 while the test profiles, processed in an 
identical manner, did not. This appears to 
rule out the possibility that our reduction 
process was responsible for the line narrow- 
ing above n = 200 since it did not affect 
the test lines. Further proof of this claim 
is presented in §4 below. 

As was found for the observations be- 
low ~ 3 GHz, the 6 GHz data also show 
an unexplainable line narrowing. This 



mysterious result has since been exam- 



ined by several au t hors (iHev 12006 



2009; Watson 2006; Oks 2004; Griem 



Gigosos et al. 




Gordon 



2007 



2008 



2007 . 



2005 



Gavrilenko and Oki 



von Prochazka et al.l 



but all attempts to explain it as an 
impact narrowing for quantum numbers 
above n = 200 and above An = 8 have 
been unsuccessful. Although the reduced 
impact broadening below 3 GHz might 
be explainable by variable density mod- 
els, this is unlikely to be the explanation 
for the decrease seen at 6 GHz since these 
observations were all made with the same 
telescope, the same beamwidth, and near 
the same frequency. 

2.2. Data at 17.6 GHz 

Because of the correlation between n 
and An, it is then also true that the line 
narrowing observed may be tied to An in- 
stead of n since a similar result will be ob- 
tained if the linewidths in Fig 1 are plotted 
versus An. In this case the line narrowing 
will begin at An = 8. However, in Fig 2 the 
linewidths observed for Orion at 17.6 GHz 
are plotted versus n. The numbers give 
the An-value associated with each transi- 
tion. Since the n-values do not go above 
200 they give no new information. How- 
ever, here the widths of the lines show no 
sign of narrowing even though An-values 
as high as An = 17 have been detected. 
This means that we cannot argue that the 
line narrowing is tied to a given An-value, 
an d this agrees w ith what was concluded 
by IWatsonI (120061 ). Because of this, it also 
seems unlikely that the narrowing is tied 
to a particular n-value. What is required 
is a parameter that can be shown to fit all 
the data, explaining not only 1), the line 
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narrowing above (202,8) at 6 GHz, but also 
2), the lack of line narrowing for An lines 
above 8 at 17.6 GHz, as well as 3), the nar- 
rowing at observing frequencies below ~ 3 
GHz. 

3. Density of Adjacent n-Transitions 
in Frequency Space 

One of the most obvious properties of 
RRLs is the fact that the spacing be- 
tween adjacent n-transitions decreases sig- 
nificantly as the observing frequency de- 



tected by ISmirnov et al.l (119841 ) at 5 GHz, 
the open c ircles represent t he transitions 
detected bv lSell et all fl201l[ ) at 6 GHz, the 



filled squares those detected by iBell et al. 



from recombination line tables 


Lillev and Palmer 


1968; 


Towle 


1993; 


Towle et al. 


1996). Here 



we define the density of the n-transitions, 
Dni in number of lines per GHz, as 

Dn = 1000/[z/(„^An) - l^{n+lAn)\ 

where v is the frequency of the transi- 
tion in MHz. 

It should be noted that Dn is much 
smaller than Dt-, where Dt is the density 
of all lines in a given frequency interval. 
Unlike Dt does not vary with An for 
a fixed observing window. 

In Figs 3 and 4, D^ is plotted versus 
n and An respectively for several different 
observing frequencies indicated by the la- 
bels. First, it is clear from these plots that 
the line density, D„, varies significantly 
with both n and An, as well as with the 
observing frequency. The horizontal line 
at Dn = 11.6 in each figure represents the 
density at which impact broadening at 6 
GHz was observed to begin to decrease and 
corresponds to n = 202 and An = 8. 

In Fig 4, the An transitions detected 
by different investigators at several differ- 
ent observing frequencies have been indi- 
cated (the filled dots are the transitions de- 



(I2OIII ) at 17.6 GHz, the open squares are 
detections made by iRood et al.l (119841 ) at 
8.7 GHz). The stars indicate the loca- 
tions of the Hna transitions observed by 
early investigators. Only for line den- 
sities above Dn = 11.6 was the impact 
broadening less than predicted. Recently, 
the GET was used to detect Voigt pro- 
f iles in several Hnct lines iii Sagittarius B2 
(jvon Prochazka et al.ll2010l ). Since the Dn 
values of these lines were all less than 10 
they cannot give us any new information 
on line narrowing mystery. 

It is clear that if Dn = 11.6 is the density 
at which the impact effect begins to be re- 
duced, the 6 GHz data would show the line 
narrowing while the 17.7 GHz data would 
not, even though lines with An- values up 
to 17 were observed. This then can explain 
the first two mysteries. 

Finally, in Fig 3 the filled circles repre- 
sent the Una transitions observed at each 
of the different observing frequencies. The 
D^-values for these points have been re- 
plotted versus observing frequency in Fig 
5. Here it is clear that for low-An tran- 
sitions, if it is true that the impact effect 
is reduced above Dn = 11.6, a reduction 
in line width is predicted to begin below 
~3.5 GHz. This is exactly what was found 
by the early investigators and would mean 
that some process other than the proposed 
variable density model might be the real 
explanation for their observations. 

The increasing line density model can 
thus be tied to the effects observed in all 
of the above three cases. 
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4. Discussion 

Although it has been demonstrated 
above that the hne narrowing at 6 GHz 
might be an effect related to the density 
of n-transitions, other possibilities need to 
be considered. We now consider whether 
a), all detections claimed are secure detec- 
tions and the line narrowing is somehow 
related to the collision process through the 
n-transition density D„, or whether b), all 
detections claimed are secure detections 
but the narrowing is somehow related to 
the data processing technique, or whether 
c), the high- An lines that are narrower 
than predicted are spurious detections. 

Previously we have concluded from the 
way the line areas in Orion and W51 arc 
observed to fall off smoothly and as ex- 
pected with increasing An, that the high- 
An transitions are valid detections, which 
immediately rules out possibility c). The 
questions that remain to be answered are 
then: What causes the line narrowing? 
Does it occur in the source, or in the data 
reduction process? With increasing line 
density now seeming to play a role it is le- 
gitimate to ask whether this line narrowing 
could somehow be caused by the reduction 
process, possibly related to line confusion, 
and we now consider this possibility. 

First, as noted above, the total line den- 
sity Df for all An values must be much 
higher than the density of the n-transitions 
related to one An value. At a given fre- 
quency the total line density for all An 
lines between 1 and 20 is easily obtained 
by counting the lines listed in recombina- 
tion line tables. For 6 GHz this value is 
~ 263 lines per GHz, which corresponds to 
an average of 1 line every 3.8 MHz. For 



linewidths close to 0.5 MHz this would not 
normally lead to confusion. However, if 
An values greater than 20 are included, 
and the observed lines are immersed in a 
low-level background of faint lines, confu- 
sion could be a problem. A simple anal- 
ysis shows that most of the noise from 
this background is at low spatial frequen- 
cies, unlike continuum noise, which could 
selectively bury the low spectral frequen- 
cies from the Voigt profile wings. The de- 
tectable signal would thus be dominated by 
high spatial frequencies, which are primar- 
ily in the central, narrow part of the lines, 
resulting in only the "tip of the iceberg" as 
it were, being detected. 

Although Dt is constant and does not 
vary with An, the lines do get weaker with 
increasing An and the wings get wider. If 
this is the correct interpretation it would 
lead to a significant loss of power in the 
lines as An increases, as well as the ob- 
served decrease in the lincwidth. Further- 
more, this effect would be present only for 
the processed telescope lines and absent 
for the processed test hues where there are 
no confusing lines, so the line narrowing 
would also be expected to show up only in 
the telescope lines as seen in Fig 1. 

Whether or not this is the correct in- 
terpretation can easily be checked by com- 
paring how the power in the line changes 
with increasing An for both the telescope 
lines and the test lines. If the line nar- 
rowing with increasing An is due to this 
confusion-related process and only the up- 
per tip of the line is being detected, the 
power detected in the lines must then de- 
crease faster for the telescope lines than for 
the test lines. 

We have previously plotted the line area 
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versus n on logarithmic scales for both the 
Orion and W51 data and found that the 
opposite is true. The area of the telescope 
lines falls off m ore slowly than the area 
of the test lines (IBell et al.ll2000l . see their 
Figs 3 and 5). This indicates that the line 
confusion process cannot be the cause of 
the line narrowing, leaving us with the like- 
lihood then that the increasing line density 
model discussed above is the true explana- 
tion for the line narrowing. 



In our previous paper (iBell et al.ll2000[ ) 

we expressed concern over why the pro- 
cessed telescope line areas remained above 
the processed test line areas while the line 
widths narrowed for An > 8 at 6 GHz. 
There is now a simple explanation for this 
mystery. The da shed line in Figs 3 and 5 of 
Bell et al.l ( 120001 ) show how the true areas 
of the lines are expected to decrease with 
increasing n and An. The solid curve in 
Figs 3 and 5 shows how the processed line 
area is reduced for the test lines by the 
overlapping process associated with fre- 
quency switching. This reduction increases 
with An because the linewidth increases. If 
something causes the line narrowing to oc- 
cur before the data reduction takes place, 
but without loss of power in the lines, since 
narrower lines are less affected by the over- 
lapping process the line areas would then 
be predicted to fall off slower than the test 
lines, as found. This slower fall-off is then 
exactly the evidence required to confirm 
that the linewidth is truly increasing less 
rapidly with An in the source lines than 
in the test lines, which indicates that the 
line narrowing process must be source re- 
lated, taking place prior to the data reduc- 
tion process. 

It is concluded that all the evidence then 



indicates that the line narrowing seen in 
impact broadened recombination lines is 
source related occurring prior to the data 
reduction stage. The evidence also shows 
that it can be linked directly to the increase 
in the n-transition density Dn that occurs 
as the observing frequency decreases. 

Because this effect is tied to the density 
of adjacent n-transitions, as opposed to the 
total line density in frequency space, this 
may suggest that the effect is related to 
the fact that in a collision it is only those 
n-levels that are close to one another that 
will be affected while more distant n-levels 
are not involved, even though they may be 
part of transitions that produce lines in the 
same frequency space. 

Finally, we might also ask if the ef- 
fect we are discussing could be due to 
masering. This would appear to be un- 
likely. The level populations and their in- 
versions at high n are well studied and ap- 
pear to rule out strong maser amplifica- 



tion at cm wavelengths f lStrelnitski et al. 



19961 ). Strong line masering should occur 



only at mm wavelengths in regions of large 
mm optical depth where the combined line 
and continuum optical depth can be large 
and negative. These conditions are rare 
and are convincingly seen in only one or 
two sources (e.g. MWC 349). The lines 
are also variable. At cm wavelengths, the 
net optical depths (line -|- continuum) are 
always positive, preventing strong maser 
lines. There can be a modest enhance- 
ment of the line intensity which is pro- 
duced by a reduced absorption of the asso- 
ciated continuum at the line centre, and 
which then appears as a slight enhance- 
ment in net emission at the line centre over 
that of LTE. The enhancements are usually 
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10 — 30% but there is no line narrowing. 



5. Conclusion 

The instances where impact broadened 
linewidths have been found to be narrower 
than expected have been re-examined. It 
is demonstrated that the line narrowing for 
high An values at 6 GHz must have oc- 
curred prior to the data processing stage 
and therefore must be source related. It 
is also demonstrated that it is unlikely to 
be due to line confusion. It is further 
demonstrated that the observed line nar- 
rowing can be tied to one of the funda- 
mental properties of radio recombination 
lines; namely the fact that the spacing of 
adjacent n-transitions increases with fre- 
quency. The line narrowing begins when 
the n-transition density, D„, exceeds ~ 
11.6 transitions per GHz. This may im- 
ply that it is somehow related either to a 
previously overlooked effect in the impact 
broadening process, or to some unknown 
parallel process, that is also tied to the n- 
transition line density. Based on these re- 
sults it can also be concluded, in agreement 
with several previous theoretical investiga- 
tions, that the line-narrowing effects are 
not tied to any particular range of either 
n or An. 
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Fig. 
from 



1. — Obser v ed lin e width at 6 GHz 
Bell et al.l feoilh after 



( 1201 ih alter processing 
plotted vs n for Orion (filled circles) and 
W51 (open squares). The numbers rep- 
resent the relevant An values. The curve 
represents the line widths obtained using 
the Voigt profile test lines generated us- 
ing impact broadening theory and known 
source parameters. There is no evidence 
that the reduction process has affected the 
width of the test hues above n = 200. Here 
n is the value of the lower level of the tran- 
sition. An = 8 corresponds to n = 202 in 
this figure because (n, An) = (202,8) is the 
An = 8 transition whose frequency falls in- 
side our observing window. 



Fig. 2. — Observed line width at 17.6 



GHz from iBell et al.l ( 1201 ll ) after process- 
ing, plotted vs n for Orion. Numbers are 
the associated An-values. They show no 
sign of narrowing above An = 8 which is 
indicated by the dashed line. 
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Fig. 3. — Plot showing how changes 
with n for different observing frequencies. 
Impact broadening is reduced above a hne 
density of 11.6 hues per GHz which occurs 
at n = 200 for hues near 6 GHz. Thus no 
narrowing would be expected at 17 GHz, 
while little broadening would be expected 
at 1.4 GHz. The filled dots represent Hna 
lines at each frequency. 




Fig. 4. — Plot showing how Dn changes 
with An. (filled circles) the transitions 
detected at 5 GHz from Smirnov et al. 



(jl984j ). (ope n circles) t he 6 GHz tran 



sitions from iBell et all (1201 If ). (filled 



squares) the 17 .6 GHz transitions from 
Bell et al.l (l201l[ ). (open squa res) the 8.7 



GHz transitions detected by iRood et al. 



(Il984j ). (stars) Hna lines. If impact broad- 
ening begins to be reduced above D„ = 
11.6 it would only be apparent so far in 
the 6 GHz data and the low frequency Hnct 
hues. 
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Fig. 5. — Plot showing how the hne den- 
sity for alpha lines changes with observing 
frequency. Impact broadening is reduced 
above D„ = 11.6 lines per GHz which oc- 
curs for Hydrogen alpha lines observed be- 
low ~ 3.5 GHz, as found by early investi- 
gators. 
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